Yoctosecond metrology through HBT correlations from a quark-gluon plasma 
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Expansion dynamics at the yoctosecond timescale aflFect the evolution of the quark gluon plasma 
(QGP) created in heavy ion collisions. We show how these dynamics are accessible through Hanbury 
Brown and Twiss (HBT) intensity interferometry of direct photons emitted from the interior of the 
QGP. A detector placed close to the beam axis is particularly sensitive to early polar momentum 
anisotropics of the QGP. Observing a modification of the HBT signal at the proposed FoCal detector 
of the LHC ALICE experiment would allow to measure the isotropization time of the plasma and 
could provide first experimental evidence for photon double pulses at the yoctosecond timescale. 
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The quark-gluon plasma (QGP) as created in heavy 
ion colliders like RHIC or LHC exists for a duration of 
a few tens of yoctoseconds (1 ys = 10~^^s). During its 
expansion, the plasma is exposed to different momen- 
tum anisotropies. Azimuthal anisotropies arise in non- 
central collisions and are responsible for the elliptic flow 
[l|, Polar anisotropies arise at very early times right 
after the collision due to the longitudinal expansion of 
the plasma 0|. The latter kind of anisotropies causes 
a variety of fascinating effects: early polar momentum 
space anisotropies can induce Chromo-Weibel plasma in- 
stabilities y48| , could allow for a violation of the viscosity 
bound 0, or can lead to photon double pulses that are 
separated merely by yoctoseconds [lo| . 

Photons are a particularly suitable probe for the early 
phase of the plasma, because once they are produced 
through quark Compton scattering or quark-antiquark 
annihilation, they leave the strongly interacting plasma 
likely without further interaction. It turns out that the 
photon production process is strongly polarization and 
direction dependent [ill . |12| . A strong polar anisotropy 



can lead to temporary suppression of photon emission in 
forward direction and thus to non-trivial pulse shapes 
that differ from the decay one would expect from an 
isotropically cooling plasma. Under particular condi- 
tions, even double pulses seem possible Measur- 
ing the pulse envelope of photons on the yoctosecond 
timescale would therefore provide firsthand information 
about this early evolution. 

Unfortunately, there are no detectors available yet that 
could time-resolve a possible signal at the yoctosecond 
scale. State-of-the-art laser physics deals with attosec- 
ond metrology [l3|, . The next generation of laser fa- 
cilities like the Extreme Light Infrastructure (ELI) flB l 
or the IZEST [3 

strive to produce zeptosecond photon 
pulses. Even though there are suggestions to characterize 
photon pulses down to zeptosecond timescale [l7|, their 



applicability to photons from the QGP seems doubtful. 
A feasible way to resolve the space-time dynamics at the 
femtometer and yoctosecond scale is the Hanbury Brown 



and Twiss (HBT) interferometry. Different from intensity 
interferometry of hadrons (like tt or 77) which essentially 
probe the surface of the plasma, photons provide infor- 
mation from the interior of the plasma [l8l.ll9j. Previous 



calculations of photon interferometry for central |2d . |21 1 
and non-central [22] heavy- ion collisions, including effects 
of extremely strong magnetic fields [23|, have assumed 
isotropic photon emission even at the very early stages 
of the QGP. 

In this Letter, for the first time, we take into account 
early polar momentum space anisotropies for the pho- 
ton emission rates to calculate two-photon momentum 
correlations. Modifications of the correlation function of 
photons emitted close to the beam axis can be linked to a 
non-trivial temporal evolution of the photon emission due 
to polar anisotropy. We show that the detection of inten- 
sity correlations close to the beam axis, as illustrated in 
Fig. [ll will be feasible in the future Forward Calorimeter 
(FoCal) detector which is likely to be installed in 2017/18 



during the ALICE detector upgrade |2J|. Our method 
would allow to measure the isotropization time of a QGP, 
and could also establish experimentally the existence of 
photon double pulses at the yoctosecond timescale. 

The photons of a few GeV energy we want to detect 
are predominantly produced at very early times, within a 
few yoctoseconds after the collision. In the description of 
correlation functions, we can thus neglect contributions 
from radiative decays of long-lived hadrons like 7r° and 
77, because of much larger length and time scales involved 
which translate to relative momenta that are smaller by 
orders of magnitude [2^. Also, we neglect a transverse 
expansion of the system, since the high-energy photons 
are most likely emitted close to the center of the QGP. 

The HBT correlation function for two photons 
with momenta k and k' is given by C2(k, k') = 
P2(k,k')/Pi(k)Pi(k'), where Pi(k) = Jd'^xwix.k) is 
the single-particle and P2IIL '^') the two- particle inclu- 
sive distribution function 12511 
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The factor i is a sta- 



tistical factor from averaging over the photon spin |26| . 
The source function ■w{x,k) = dR{x,'k)/ {d'^x d^k) de- 
scribes the mean number of particles of four- momentum k 
emitted from a source element centered at the space-time 
point X [Hi. a chaotic source is assumed, as contribu- 
tions by correlated two-photon emissions are estimated 
to be negligible (20| . 

In the framework of the current formalism (l8| one can 
show analytically that a source function composed of two 
temporally separated Gaussians leads to oscillations in 
momentum space in the HBT functions. However, oscil- 
lations in the correlation functions could also be caused 
by two spatially separated emission centers. The ques- 
tion arises how to distinguish a temporal variation from 
a spatial variation in the HBT signal. This is possible 
by combining correlation measurements in different di- 
rections: HBT oscillations due to temporal separation 
are independent of the direction of observation, while 
spatially separated sources lead to a strong directional 
dependence of the HBT signal. 

In practice, emission centers will not follow a perfect 
Gaussian shape. It has been found that the two-photon 
correlator shows a strongly non-Gaussian shape along the 
polar direction [2^ . In the azimuthal direction, oscilla- 
tions can appear due to a non-circular intersection region 



in non-central collisions (22| . In the following we present 



evidence that early polar momentum space anisotropics 
can lead to observable modifications of the HBT sig- 
nal, including the appearance of a side peak, although 
we are aware that it will be experimentally challeng- 
ing to distinguish such modifications from other possi- 
ble sources of oscillations in the HBT correlation func- 
tions. Note that there is also the possibility of introduc- 
ing fake oscillations due to inaccurate numerical integra- 
tions 2^, 23 ■ We therefore carefully cross-checked our 
numerical results presented below using analytical and 
semi-analytical models of temporally separated sources. 

The longitudinal expansion of a plasma right after a 
collision can be described in two limiting cases: one 
is the ideal hydrodynamical evolution as described in 
the Bjorken expansion picture [2^, where quark and 
gluon distribution functions stay isotropic throughout 
the expansion. The other extreme is the free stream- 
ing limit (29I . which neglects all parton interactions, 
and where momentum anisotropy increases over time. 
We base our calculation on a model that can interpo- 
late between these two limiting cases. The momen- 
tum anisotropy is implemented through a modification 
of isotropic distribution functions /iso according to [13 



/(p) = /iso(v/p' + e(p-n)2), (2) 
where n points along the beam axis and the anisotropy 




Figure 1: (color online) Concept of HBT detection. Two 
heavy ions collide with impact parameter b and produce a 
QGP that rapidly expands at early times along the beam 
axis (z-axis). Photons are emitted from the QGP and arrive 
at a detector placed in forward direction, like the proposed 
FoCal detector [2^ . HBT correlations are obtained from pho- 
tons that arrive simultaneously, whose momentum vectors are 
separated by an angle 0^,,i. 



parameter ^ = (p|,)/(2(p|)) — 1 is defined in the range 
— 1 < ^ < 00. Values of ^ > contract an isotropic distri- 
bution along the beam axis, while ^ < stretch it. A sin- 
gle parameter S can describe the scaling solutions of ideal 
hydrodynamical evolution {S — 0) and free-streaming 
expansion (S = 2), as well as other expansion scenar- 
ios like momentum-space broadening due to interactions 
((5 = 2/3) [U. The time evolution of the anisotropy pa- 
rameter ^ or the hard momentum scale Phard, which plays 
the role of temperature in an anisotropic medium, is then 
given by 



1, 



Phard(T) = — 



(3) 
(4) 



We follow the model of the plasma evolution of Ref. (32| 
who introduced a smeared step function A(r, riso,7). By 
basically replacing 5 5(1 — A(r, riso, 7)), this function 
governs the change of 5 from e.g. 2 or 2/3 to (5 = at 
approximately the isotropization time Tjso with a param- 
eter 7 which determines the smoothness of the transition. 
The initial temperature distribution for central and non- 
central collisions is assumed to be proportional to the 
thickness functions of the colliding nuclei according to 
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Figure 2: (color online) HBT photon correlation for a coUinear 
configuration. The detector is placed at a fixed angle 6 — 25° 
{rj ~ 1.5) and the difference q = k' — k between coUinear 
photon momenta k' || k with k = 4GeV/c {kr ~ 1.7GeV/c) 
is varied. The various isotropization timescales correspond 
to initial free streaming expansion with intermediate polar 
momentum anisotropy starting from nso = 2fm/c (solid line) 
down to Tiso = To (dotted line) which corresponds to an ideal 
hydrodynamic expansion. The impact parameter is chosen as 
6 = 10 fm and the anisotropy model parameter as 7 = 2. 



Figure 3: (color online) HBT photon correlation for a detector 
placed in forward direction. One momentum vector is placed 
at 2° away from the beam axis {ri « 4) and the other varied 
between 2° and 8° (r; « 4 to « 2.7, at the same azimuthal 
angle) at fixed momentum k — 25GeV/c. A distinct differ- 
ence and change of shape of the HBT function with respect to 
the isotropization time of the quark-gluon-plasma is observed. 
The impact parameter is fe = 12 fm and the anisotropy model 
parameter is 7 = 2. 



the Glauber model [33|. As in Ref. [3J|, we use hard 
spheres to model the initial nuclear charge density. 

The main contributions of photon production arise 
from quark-Conipton scattering and quark-antiquark an- 
nihilation processes [sBj. In an anisotropic plasma, the 
corresponding photon production rate EcPR/d^k can 
only be obtained numerically and shows strong direc- 
tional dependence [III 12 1- Corresponding expressions 



for bremsstrahlung or inelastic pair annihilation become 
important at lower energies. Therefore, as in Ref. [ll|, 
we do not take into account these soft scattering pro- 
cesses. As in Ref. [l3|, we use parameters that are rel- 
evant to heavy ion cohisions at the LHC, with initial 
temperature Tq = 845 MeV, plasma freezeout tempera- 
ture Tc = 160 MeV, nucleus radius R = 7.1 fm, and a 
plasma formation time of tq = 0.088 fm/c. To estimate 
the effect of early polar momentum anisotropics, various 
isotropization times starting from ideal hydrodynamic 
expansion (tiso = tq) up to ideal free-streaming 6 = 2 
with isotropization time Tjso = 2 fm/c are compared. 

We consider two different configurations for the calcu- 
lation of the HBT correlation function: a coUinear and 
a non-collinear configuration. In the coUinear configura- 
tion in Fig. [21 the two photon momentum vectors k and 
k' point into the same direction (k [| k') at a polar angle 9 
away from the beam axis. As in |l9l.l25|. the free param- 
eter of the two-particle HBT function is the momentum 



difference q = k' — k. We focus on collisions with a large 
impact parameter b > 10 fm so that a possible signal is 
not averaged out by the transverse size of the plasma [l3| . 
For hydrodynamic expansion (riso = Tq), the correlation 
function shows the unobtrusive behavior of a monoton- 
ically decreasing function. If one includes the effect of 
early polar momentum anisotropics with Tiso ^ 1.5 fm/c 
however, the HBT function reveals a non-trivial shape 
and exhibits a plateau-like structure. The reason for the 
appearance of such plateaus is the modification of the 
photon production rate in an anisotropic plasma. The 
suppression of the photon production in forward direc- 
tion at times tq < t < Tjso before isotropization leads to 
a non-trivial emission envelope which could result in two 
temporally separated peaks jlO| . Such an emission enve- 
lope leads to a side peak in the correlation function. A 
drawback of this configuration is that coUinear photons 
can not be readily resolved by photon calorimeters. 

Therefore we also consider a non-collinear configura- 
tion as shown in Fig. [31 In this case, both photon mo- 
menta share the same magnitude (|k| = |k'|), but are 
positioned at a relative angle 0roi to each other within 
the reaction plane. The selected parameter range is cov- 
ered by the proposed FoCal detector In this forward 
detector, photons arrive highly blue-shifted. A fixed pho- 
ton momentum fc = 25 GeV/c as observed by the detector 
corresponds to a transverse momentum which decreases 
from kr ~ 3.5GeV/c at ry « 2.7 to kr « 0.9GcV/c 
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at ry ~ 4. Thus, these photons are most hkely emit- 
ted from the early QGP. We see similar behavior in the 
non-collinear HBT function as in the collinear configura- 
tion. Early polar momentum space anisotropics result in 
a narrower correlation function as well as the emergence 
of a second maximum. In principle, two temporally sep- 
arated photon emission peaks would lead to an oscilla- 
tion of the correlation function, but beyond the first two 
peaks, further maxima are not identifiable for LHC pa- 
rameters after integration over the space-time evolution 
of the QGP. The distance between the main peak and the 
side peak of the correlation function is inversely propor- 
tional to the time interval between two temporally sepa- 
rated peaks in configuration space, and thus also roughly 
inversely proportional to the isotropization time. An as- 
sumed isotropization time of riso = 2 fm/c leads to peaks 
in the correlation function separated by about 3° for the 
parameters chosen. Such a structure can in principle be 
observed in the proposed FoCal detector which will be 
able to resolve photons that are separated merely by a 
fraction of a degree [2J|. 

The photon momentum correlations presented here 
could be influenced by various effects: Although trans- 
verse expansion can be neglected at very early times close 
to the center of a collision, the question is more delicate 
for highly non-central collisions. The space-time evolu- 
tion in transverse direction may produce additional mod- 
ifications of the photon correlation functions. Also, at 
very large rapidity one can not assume a boost-invariant 
particle multiplicity. Taking this effect into account will 
affect the absolute rate of the observed photons, but it 
will not destroy the correlation between them. 

Regarding the feasibility of the detection, it will be 
experimentally challenging, but not impossible. If one 
assumes an annual yield of at least 10^ prompt photons 
in the energy range kx = I GeV/c to 4 GeV/c from heavy 
ion collisions that hit the FoCal detector [36|, one would 
observe a few hundred photon pairs within the same time 
frame. The modification of the signal presented here is 
only caused by direct photons and should therefore be 
distinguishable from background photons. 

To summarize, we have calculated the intensity cor- 
relation of photons produced at an early stage of the 
quark-gluon plasma, taking into account full polar mo- 
mentum anisotropy. Besides other known sources that 
could lead to oscillations in the HBT signal, we found 
that the correlation function is particularly sensitive to 
the early time evolution of the plasma. Large isotropiza- 
tion times lead to distinctive modifications of the HBT 
correlation function with a side peak appearing a few 
degrees separated from the main peak. The detection 
of such a structure, for example in the proposed FoCal 
detector at the ALICE experiment, would enhance our 
knowledge about the early evolution of the QGP, includ- 
ing information about the isotropization process and the 
isotropization time. It could also provide first indirect 



experimental evidence for possible photon double pulses 
at the yoctosecond timescale. 
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